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Abstract 

We consider a general N-degree-of-freedom dissipative system which admits 
of chaotic behaviour. Based on a Fokker-Planck description associated with 
the dynamics we establish that the drift and the diffusion coefficients can be 
related through a set of stochastic parameters which characterize the steady 
state of the dynamical system in a way similar to fluctuation-dissipation re- 
lation in non-equilibrium statistical mechanics. The proposed relationship is 
verified by numerical experiments on a driven double well system. 
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I. INTRODUCTION 



Although deterministic in principle classically chaotic motion is stochastic in nature. 
Ever since the early numerical study of Chirikov mapping |]J revealed that the motion of a 
phase space variable can be characterized by a simple random walk diffusion equation, at- 
tempts have been made to describe the chaotic motion in terms of Langevin or Fokker-Planck 
equations PIH. It is therefore easy to comprehend a close connection between classical chaos 
and statistical mechanics. Two distinct situations arise in this context. The first one con- 
cerns whether classical chaos may serve as a basis for classical statistical mechanics since the 
ultimate justification of the postulates of statistical mechanics like Boltzmann hypothesis 
of molecular chaos, ergodicity or the postulate of equal a priori probability rests on the 
dynamics of each particle [0-0]. The second one concerns the following : Given that the 
classical chaotic motion is stochastic how and to what extent one can realize the formulation 
of statistical mechanics for useful description of classical chaos P-f2~H| keeping in mind that 
one essentially deals here with a few-degree-of-freedom system. The present paper addresses 
the second issue. 

The emergence of stochastic behaviour of the classically chaotic system is due to the loss 
of correlation of initially nearby trajectories. This is reflected in the nature of the largest 
Lyapunov exponent [32] whose calculation rests on the linear equation of motion for the sep- 
aration of these trajectories. When chaos has fully set in, the time dependence of the linear 
stability matrix or Jacobian of the system [2j| in the equation of motion in the tangent space 
can be described as a stochastic process since the phase space variables behave as stochastic 



variables. In a number of recent studies we have shown [I7|-pi[ that this fluctuation of the 
Jacobian is amenable to a theoretical description in terms of the theory of multiplicative 
noise. This allows us to realize a number of important results of nonequilibrium statis- 



tical mechanics, like Kubo relation ||17|| , fluctuation-decoherence relation ||18|| , exponential 
divergence of quantum fluctuations |T9|-f21|, thermodynamically inspired quantities, e. g. , 
entropy production in chaotic dynamics. Based on a Fokker-Planck description in the tan- 
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gent space where the drift and the diffusion coefficients explicitly depend on the phase space 
variables or dynamical properties of the system, we show that a connection between the two 
moments in terms of the stochastic parameters which characterize the long time limit of the 
dynamical system can be established in the spirit of fluctuation-dissipation relation. We 
verify the theoretical proposition by numerical experiments on a simple dissipative system. 

The rest of the paper is organized as follows : In Sec. II we introduce a Fokker-Planck 
description of the dynamical system in the tangent space and identify the drift and diffusion 
coefficients which are the functions of fluctuations of the phase space variables. This is fol- 
lowed by solving the Fokker-Planck equation for the steady state distribution required for the 
calculation of long time averages in Sec III. In Sec. IV the dynamical stochastic parameters 
which characterize the long time behaviour of the system are introduced. The first one of 
them is a well-known stochastic parameter closely related to Kolmogorov entropy. With the 
help of these stochastic parameters we establish a connection between the drift and diffusion 
coefficients of the Fokker-Planck equation in the spirit of fluctuation-dissipation relation in 
nonequilibrium statistical mechanics. In Sec. V we illustrate the general method by an 
explicit numerical example to verify the theoretical proposition. The paper is concluded in 
Sec. VI. 

II. A FOKKER-PLANCK EQUATION FOR DISSIPATIVE CHAOTIC 

DYNAMICS 

We are concerned here with a general N-degree-of-freedom system whose Hamiltonian is 
given by 

N 2 

H = T,£r + v (tiiU) , i = i~.N (l) 

i=l % 

where {q%-,Vi\ are the co-ordinate and momentum of the i-th degree-of-freedom, respectively, 
which satisfy the generic form of equations 

. OH , . dH 

dpi dqi 
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We now make the Hamiltonian system dissipative by introducing — jpi on the right hand 
side of the second of Eqs.(2). For simplicity we assume 7 to be the same for all the N degrees 
of freedom. By invoking the symplectic structure of the Hamiltonian dynamics as 

qi for i = 1 • • • N , 
p^ N for % = N+ 1,---2N . 

and defining I as 

E 
-E -7E 

where E is an TV x N unit matrix, and is an N x N null matrix, the equation of motion 
for the dissipative system can be written as 

z * = L%a7- ( 3 ) 
3=1 oz i 

We now consider two nearby trajectories, Z{, z% and Zi + Xi, Z{ + Xi at the same time t in 
2N dimensional phase space. The time evolution of separation of these trajectories is then 
determined by 



2N 



(4) 



in the tangent space {Xi}, where 



Jij — I; 



ik 



d 2 H 

dzudz- 



(5) 



Therefore the 2N x 2N linear stability matrix J assumes the following form 



J 



E 
M(t) -^E 



(6) 



where M is an iV x matrix. Note that the time dependence of stability matrix J_(t) is 
due to the second derivative q^q z . which is determined |^| by the equation of motion (3). 
The procedure for calculation of and the related quantities is to solve the trajectory 
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equation (3) simultaneously with Eq.(4). Thus when the dissipative system described by 
(3) is chaotic, J(t) becomes (deterministically) stochastic due to the fact that Zi-s behave as 
stochastic variables and the equation of motion (|J) in the tangent space can be interpreted 
as a stochastic equation fi~7^pTf. 



In the next step we shall be concerned with a stochastic description of J(i) or M(t). For 
convenience we split up M into two parts as 

M = Mp + M 1 (t) (7) 

where Mq is independent of variables {z{\ and therefore behaves a sure or constant part and 
Mi is determined by the variables {z;} for i = 1 • • -2N. Mi refers to the fluctuating part. 
We now rewrite the equation of motion (f|) in tangent space as 

X = JX 

= L({X i },{z i }) (8) 

where X and L are the vectors with 2N components. Corresponding to (0) L in can be 
split up again to yield 

X = L°(X)+L 1 (X,{z i (t)}) , i = l---2N . (9) 

Eq.(|4]) indicates that Eq.([8]) is linear in {Xi}. Eqs.(|4]), (||) and (||) express the fact the 
first iV components of L 1 are zero and the last iV components of L 1 are the functions of 
{Xi} for i = 1 • • ■ N. The fluctuation in L\ is caused by the chaotic variables {^}-s. This 
allows us to write the following relation (which will be used later on), 

V x • £V(TO) = L 1 ■ V x 0({X,}) (10) 

where 4>({Xi}) is any function of {^}. Vj refers to differentiation with respect to compo- 
nents {X^ (explicitly X { = Ag< for i = 1 - • • N and Xi = Aft for % = N + 1 • ■ • 2N). 

Note that Eq.(^) by virtue of (§) is a linear stochastic differential equation with multi- 
plicative noise where the noise is due to {z,i} determined by equation of motion (H). This is 
the starting point of our further analysis. 



Eq.(P) determines a stochastic process with some given initial conditions {Xj(0)}. We 
now consider the motion of a representative point X in 2N dimensional tangent space 
(X\ ■ ■ ■ X 2 n) as governed by Eq.(|9|). The equation of continuity which expresses the con- 
servation of points determines the variation of density function <p(X, t) in time as given 
by 

d<j>{X,t) 



dt 

Expressing A and A x as 



-V*-L(t)4>{X,t) . (11) 



A = -V x ■ L° and A x = -V x • L 1 (12) 
we may rewrite the equation of continuity as 

?^^ = [A + aA 1 (tMX,t) . (13) 

It is easy to recognize that while A Q denotes the sure part Ai contains the multiplicative 
fluctuations through {zi(t)}. a is a parameter introduced from outside to keep track of the 
order of fluctuations in the calculations. At the end we put a = 1. 

One of the main results for the linear equations of the form with multiplicative noise 



may now be in order [25]. The average equation of (0) obeys [ P(x,t) 



P = {A + a(A 1 ) + a 2 J™ drUA^t) exp(ry4 )A 1 (t - r))) exp(-rA )| P{x, t) . (14) 

The above result is based on second order cumulant expansion and is valid when fluctu- 
ations are small but rapid and the correlation time r c is short but finite or more precisely 

((A 1 (t)A 1 (f)))=0 for \t-t'\ > r c (15) 

We have, in general, (A x ) ^ 0. Here ((• • •)) implies ((&<?)> = (dCj) ~ (G){(j)- 

The Eq . (|T^) is exact in the limit r c — > 0. Making use of relation (0) in (p!4]) we obtain 



dP 
~dt 



|-V ■ L° - a(V ■ L 1 ) + a 2 J™ dr((V ■ L l (t) exp(-rV ■ L°) 
V-L 1 (t-r)))exp(rV-L°)}P . (16) 



6 



The above equation can be transformed into the following Fokker-Planck equation (a = 1) 
for probability density function P(X,t), (the details are given in the Appendix A); 



dP(X, t) 

Ft 



(17) 



where, 



F = L° + (L 1 ) + Q 



(18) 



and Q is a 2A r -dimensional vector whose components are defined by 



Qi = -J «i^»drdi(r)tZ 2 (r) 



Here the determinants Deti(r), Det 2 (r) and i?' are given by 



(19) 



Detxir) 



dX~ T 



dX 



and Det2{r) 




dX 



dX- 



and i5 = E4(^*)^E4(-^ T ,*-T 



* k 



dXj 



(20) 



roc rj V . 

Vii= / T.^L\{X,t)Ll{X-\t-r)-^))Det 1 {r)Det 2 {r)dr 
Jo , aAj, 



(21) 



It is easy to recognize F as an evolution operator. Because of the dissipative perturbation 
we note that div F < 0. 

The diffusion coefficient in Eq. ([T7|) is defined as 

dXj 

k """fc 
We have followed closely van Kampen's approach p5 to generalized Fokker-Planck equation 
(p|). Before concluding this section several critical remarks regarding this derivation need 
attention: 

First, the stochastic process Mi(t) determined by {zi} is obtained exactly by solving 
equations of motion (|3p for the chaotic motion of the system. It is therefore necessary to 
emphasize that we have not assumed any special property of noise, such as, Mi (t) is Gaussian 
or ^-correlated. We reiterate Van Kampen's emphasis in this approach. 

Second, the only assumption made about the noise is that its correlation time r c is short 
but finite compared to the coarse-grained timescale over which the average quantities evolve. 



Third, we take care of fluctuations upto second order which implies that the deterministic 
noise is not too strong. 

Eq.flTTD is the required Fokker-Planck equation in the tangent space {X{\. Note that 
the drift and diffusion coefficients are determined by the phase space {zj} properties of the 
chaotic system and directly depend on the correlation functions of the fluctuations of the 
second derivatives of the Hamiltonian (H). 

III. THE STEADY STATE DISTRIBUTION AND THE CALCULATION OF 

AVERAGES 

In what follows we shall be concerned with the long time limit of the dynamical system. 
Thus the steady state distribution of the tangent space co-ordinates Xi(i = 1 • • -2N) are 
specially relevant for the present purpose. To make all these co-ordinates dimensionless we 
use the following transformations in Eq.(17) 

t' = u't , 



Ui = —r for i — 1 • • • N 



X ' for i = N + 1 • • • 2N , (22) 



to'd. 



where to' is a scaling constant having dimension of reciprocal of time (a possible choice is 
the linearized frequency of the dynamical system) and r' becomes a dimensionless variable. 
do is a constant (to be specified later) having the dimension of length. The resulting Fokker- 
Planck Eq.(17) reduces to 

^- vj,( ' )f+ r" ( 4' <23) 

Note that Eq.(23) is independent of do since F(X) is linear in {X{\ and T>(X) is quadratic 
in {^j}. Next we consider the stationary state of the system (|^ = 0) and make use of the 
following linear transformation ( with a2N = 1 ) 



2N 



(24) 



i=i 



in Eq.(23) to obtain the equation for steady state probability distribution P S (U) : 



_d_ 

dU 



\UP S (U)+V. 



d 2 P s 







(25) 



s 



dU 2 



ctj-s (i — 1 • • • 2^ — 1) are the constants to be determined. Here 



(26) 



and 



V S = Y^ V'ijOLiOL 



(27) 



and we disregard the time dependence of V under weak noise approximation, to treat T>' 
as a constant in the usual way. 

Putting (24) in (26) and comparing the coefficients of on both sides we obtain 2N 
algebraic equation (for ctj olzn-i and A). The set {a{\ and A are therefore known. 

The exact steady state solution, P s has the well known Gaussian form which is given by 



where iV is the normalization constant. Eq.(28) expresses the probability distribution of 
tangent space co-ordinates of the dynamical system in the long time limit. The important 
relevant quantity which measures the separation of initially nearby trajectories when the 
system has attained the stationary state can be computed by calculating the average of 
Yh=i Hi- Making use of the distribution (28) we obtain 



Note that the average as calculated above is a function of D s , A and cvs which are dependent 
on the phase space properties of the dynamical system. 




(28) 




(29) 
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IV. STOCHASTIC PARAMETERS, CONNECTION BETWEEN V s AND A ; 
FLUCTUATION-DISSIPATION RELATION 



Eq.(25) is a steady state Fokker-Planck equation in tangent space with linear drift and 
constant diffusion coefficients where the co-ordinates have been expressed as dimensionless 
variables {yi}- A and D s are the first and second moments, respectively, of the underly- 
ing stochastic process. Our objective here is seek for a connection between the two mo- 
ments. In standard nonequilibrium statistical mechanics this connection is expressed by the 
fluctuation-dissipation relation through temperature, an equilibrium parameter character- 
izing the equilibrium state. Our approach here is to follow a somewhat similar procedure. 
This implies that we search for the stochastic parameters which characterize the long time 
limit of the nonlinear dynamical system. We show that an appropriate relation between T> s 
and A can be established through these parameters. 

An important parameter proposed many years ago by Casartelli et. al. p4[ (a precursor 
for the largest Lyapunov exponent used as a measure of regularity or chaoticity of a nonlinear 
dynamical system) is the long time average of In ^ where do is the separation of the two 
initially nearby trajectories and d(t) is the corresponding separation at some time t. To 
express d{t) (having dimension of length) we write d{t) = [J2iLi(Xi) 2 + Y^=n+i(^/) 2 ]^ ■ d(t) 
is determined by solving numerically Eqs. (3) and (4) simultaneously or their appropriately 
transformed version for the initial condition zq corresponding to Eq.(3). In going from j-th 
to j+l-th step of iteration in course of time evolution any of the components of X say X{ 
has to be initialized as X\ = -f-do. This initialization implies that at each step iteration 

starts with same magnitude of do but the direction of do for step j+1 is that of d(t) for j-th 

x-' 

step (considered in terms of the ratio -^-). For a pictorial illustration we refer to Fig.l of 
Ref. ||22|| . j-th time of iteration implies t = jT (j = 1, 2 • • • oo) and T is the characteristic 
time which corresponds to the shortest ensemble averaged period of nonlinear dynamical 



system. Thus following Casartelli et. al. |23| a stochastic parameter can be defined by the 



following time average of In j- as 
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1 - d- 

a n (t,z ,d ) = -E ln X ( 3 °) 
n . do 



It has been shown |24fl that as n — * oo a n has a definite value. For the disordered system 
it is positive and for the regular system it is zero. The difference of a n from the largest 
Lyapunov exponent is also noteworthy. Our object here is to generalize (30) by defining the 
other higher order moments (higher than the first a n ->oo). To express these quantities we 
define first 

a' = ln d M (31) 
do 

We now make use of the transformation (22) to express d(t) as a dimensionless quantity 
in terms of a' as follows: 

2N 

ln5> 2 = 2a' . (32) 
i=i 

The method of cumulant expansion on the other hand tells us that the average of the 
sum of yf can be written as 

= exp f^A^ m = 1,2,3,-.. (33) 

where A m -s result from cumulants of the stochastic quantity 2a' . A m -s are calculated 
dynamically from the following relations 

A 1 =m 1 , A 2 = ^[m 2 -m 2 ] , A 3 = ^[m 3 - 3mim 2 + 2m\} 

At = ^y[ m 4 — 3ml — Amirn^ + Y2m\mi — 6mf] etc. (34) 

where nik = ^ Y%=i (^ n ^) = 1' 2> 3, 4, ...]. In the spirit of Ref. [^4| we enquire, whether 
these moments/cumulants reach their steady state values in the long time limit. We have 
numerically examined the dependence of m^-s on various parameters. The parameters are 
n, the time, do, the measure of initial separation, the characteristic time T (j-th time of 
iteration implies t = jT,j = 1, 2, • • • oo). Our observation is that the limit or limit A m 



as n — ► oo seems to exist in all cases. We have examined 0] these limits for a number 
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of test cases, e. g. , for Lorentz system, Henon-Heiles system and others. In Fig.(l) we 
exhibit a typical representative long time behaviour of the cumulants A m (m = 1 to 4) for a 
driven double well potential system discussed in the next section. It is apparent that they 
attain their long time limits as n — > oo. Secondly, the first two cumulants are much higher 
compared to others The first moment is the stochastic parameter defined by Casartelli et. 
al. |24| as a quantity closely related to Kolmogorov entropy. We are therefore led to believe 



that the quantities A m -s characterize the long time limit or the steady state of a dynamical 
system. 

The relations (33) and (29) can now be combined to give 



(35) 



The above relation is the central result of this paper. This establishes a connection 
between the drift and the diffusion coefficients of the Fokker-Planck equation (25) through 
the stochastic parameters characterizing long time behaviour of the nonlinear dynamical 
system. It must be emphasized that both the drift A and the diffusion T> s coefficients arise 
from the deterministic stochasticity implied in the dynamical equation of motion (3). The 
relation (35) is therefore reminiscent of the familiar fluctuation-dissipation relation. 

A few points regarding the relation (35) are in order. It is important to note that the 
fluctuation-dissipation relation in conventional nonequilibrium statistical mechanics is valid 
for a stochastic system for which the noise is internal. The spiritual root of this relation lies 
at the dynamic balance between the input of energy into the system from the fluctuations 
of the surrounding and the output of energy from the system due to its dissipation into 



the surrounding. The system-reservoir model p^HH developed over the last few decades 
suggests that the coupling between the system and the reservoir is responsible for a common 
origin of drift and diffusion. In the present theory this common mechanism is the fluctuations 
of the phase space variables (or second derivative of the Hamiltonian) inherent in both the 
drift A and the diffusion D s coefficients of the Fokker-Planck equation. We point out that 
the relation is still valid for the pure Hamiltonian system (7 = 0). For this reason the 
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relation (35) is somewhat formal in contrast to the standard fluctuation-dissipation relation. 



V. AN EXAMPLE AND NUMERICAL VERIFICATION 

To illustrate the theory developed above, we now choose a driven double-well oscillator 
system with Hamiltonian 



H = — + aq\ — bq 1 + eqx cos Vtt 



(36) 



where px and qx are the momentum and position variables of the system, a and b are the 
constants characterizing the potential, e includes the effect of coupling constant and the 
driving strength of the external field with frequency Q. This model has been extensively 
used in recent years for the study of chaotic dynamics flI7||T8"| , [29 . 



The dissipative equations of motion for the tangent space variables Xx and X 2 corre- 
sponding to qi and px (Eq.8) read as follows: 

r A 9l = Xx 
Apx = X 2 



d_ 

dt 



Xx 




Xx 


— J 


x 2 




x 2 



(37) 



where J as expressed in our earlier notation z\ = qx and z 2 = px is given by 



where ((t) = -Ylaz\. Eq.(| 





((t) + 2b -7 
is thus rewritten as 



V 



d_ 

di 



yX 2j 



(38) 



with 



( 



X 2 

2bXx - jX 2 



\ ( 

and L 1 = 





C(t)Xx 



where L° and L 1 are the constant and the fluctuating parts (vectors), respectively. The 
fluctuation in L 1 , i.e., in is due to stochasticity of the following chaotic dissipative 
dynamical equations of motion; 
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Z\ = z 2 and z 2 = —az\ + 2bz\ — e cos Clt — r yz 2 • (39) 

The result of Eq. ( |^5| ) can then be applied and after some algebra the Fokker-Planck equation 
(17) for the dissipative driven double- well oscillator assumes the following form: 

dP dP dP d d 2 P d 2 P 

W = -^Wl - "* Xl ex 2 + W* 2P) + v ^axMl + v *axl (40) 

where 

("OO 

V 21 = X 2 ((C(t)C(t-r)))re^dr 
Jo 

and 

V 22 = X 2 / ((C(t)C(t - r)))e-^dr - X X X 2 / «C(t)C(* " r^re^dr (41) 
Jo Jo 

with 

rod 

u 2 = 2b + c + c 2 , c 2 = / {(((t)((t-T)}}Te^ T dT and c = (0 . (42) 

Jo 

The similarity of the equation (40) to generalized Kramers' equation can not be overlooked. 
This suggests a clear interplay of chaotic diffusive motion and dissipation in the dynamics. 
Using the transformation (22) Eq.(40) can be written as 

9P dP 2 dP d . ns _ d 2 P w d 2 P 

or dyi dy 2 dy 2 oy 2 dy 1 dy£ 

where 

^ = 4 , 7 = ^7 » ^21 = ^ r((C(T / )C(T / -r)))re-^r and 

a;' 2 a/ a/ Jo 

V'22 = ^ /°°((C(r')C(r' - r)»e-^r - M%M /"((CCOC^ " r)))re^dr (44) 

and the time dependence of 7/1 and 7/2 i n the diffusion coefficients have been frozen under 
weak noise approximation. 

Now using the linear transformation (24) in Eq.(43) we obtain in the stationary state 
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(45) 



where 

U = aiyi + y 2 and XU = -a 1 y 2 - u 2 y\ + 72/2 (46) 

and 

^ s = ^'22 (47) 

where for simplicity it has been assumed that V 2 i is much small compared to the Marko- 
vian contribution V 22 . 

Comparing the coefficients of y 1 and y 2 on both sides of Eq.(46) we obtain 

A«i = — ZJ 2 and A = — a,\ + 7 

Therefore we have 



-7-V7 2 + 4u; 2 7 + ^7 2 + 4o; 2 
«! = ^ and A = . (48) 

Here the negative value of a\ is taken to make A positive for a physically allowed solution 
of the steady state distribution (49). The solution of Eq.(45) is given by 

P s = Nexp^-^-(a 1 2 y 2 + 2 ai y 1 y 2 + y 2 2 )j . (49) 

With the help of above distribution the average quantities in tangent space can be cal- 
culated. Thus we have 

^ 2 > = t(4 +1 ) • < 50 > 

The fluctuation dissipation relation (35) can then be obtained by combining (50) with 
(33) as follows ; 



5TT) exp fe A '") ' (51) 



15 



A and ot\ are to be calculated using (48). For these we require explicit numerical eval- 
uation of uj 2 as defined in (43) and (44). The dissipative chaotic motion is governed by 



equations (37) and (39). We choose the following values of the parameters [2^] a = 0.5, 
b — 10, e = 10, Q — 6.07 and 7 = 0.001. The coupling-cum-field strength e has been varied 
from set to set. We choose the initial conditions ^i(0) = —3.5 and 22(0) = which ensures 
strong global chaos. Note that C2 as expressed in (42) and in the diffusion coefficients are 
the integrals over the correlations of ((t) (C{t) is the fluctuating part of the second derivative 
of the potential V{z) and is given by ((t) = —12azf). To calculate the correlation function 
((C(t)((t — r)}) and the average (C(^)) it is necessary to determine long time series in ((t) 
by numerically solving the classical equation of motion (39). The next step is to carry out 
the averaging over the time series. For further details of the numerical procedure we refer 



to the earlier work |19|-f2T|. On the other hand the cumulants A m (m = 1, 2, 3, 4) (as defined 
in (34) and (35)) are calculated from Eqs.(37) and (39) directly. The method has already 
been outlined in Sec. (IV) and in Ref. [p^j . We then plot the theoretically calculated values 
of T> s from the evaluation of A, a\ and the cumulants for several values of the coupling con- 
stant e (Eq.36) and compare them with the diffusion coefficients obtained from the direct 
numerical integration of Eqs.(39) and (37) with the appropriate transformation (22) for the 
corresponding values of e. The result is shown in Fig. 2. It may be noted that the theoreti- 
cal and numerical results are in good agreement. The validity of the fluctuation-dissipation 
relation as proposed in Eq.(35) is therefore reasonably satisfactory. 



VI. CONCLUSIONS 

The crucial question of instability of classical motion essentially rests on the linear sta- 
bility matrix or Jacobian matrix associated with the equations of motion. While the linear 
stability analysis around the fixed points is based on the assumption of constancy of this 
matrix we take full account of the time dependence of the quantity in the chaotic regime 
by considering it to be a stochastic process, since the phase variables behave stochastically. 
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Based on a Fokker-Planck description in the tangent space we trace the origin of chaotic 
diffusion and drift in the correlation of fluctuations of the linear stability matrix. 
The main conclusions of this study are the following : 

(i) We show that a class of dynamical stochastic parameters which attain their steady 
state values in the long time limit of the dynamical system may be used to characterize 
the dynamical steady state of the system. The first one of them which was proposed by 
Casartelli et. al. |24| several years ago as a measure of the chaoticity of the system is closely 



related to Kolmogorov entropy. 

(ii) We establish a connection between the drift and the diffusion coefficients of the 
Fokker-Planck equation and the dynamical stochastic parameters in the spirit of fluctuation- 
dissipation relation. The realization of this relation in chaotic dynamics therefore carries the 
message that although comprising a few degrees of freedom a chaotic system may behave as 
a statistical mechanical system (although in a somewhat different sense). 

The theoretical relations proposed here are generic for N-degree-of-freedom chaotic 
Hamiltonian system with or without dissipation and have been verified by numerical anal- 
ysis of a driven nonlinear dissipative system. We hope that the present approach will find 
useful application in searching for the related thermodynamically inspired quantities in few- 
degree-of-freedom systems. 
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APPENDIX A: THE DERIVATION OF THE FOKKER-PLANCK EQUATION 



We first note that the operator e(— rV-L°) provides the solution of the equation [Eg. (|13|) 
a = 0] 

df(X,t) 



dt 



-V x -L»f(X,t) 



(Al) 



/ signifies the "unperturbed" part of P which can be found explicitly in terms of character- 
istic curves. The equation 



X = L°(X) 



(A2) 



determines for a fixed t a mapping from X(r = 0) to X(t), i. e., X —> X T with inverse 
(X T )- T = X . The solution of Q is 

dX-* 



f(X,t) = f(X-\0) 



dX 



e[-tV-F o ]f(X,0). 



(A3) 



d(X) 



being a Jacobian determinant. The effect of e(— tV • LP) on f(X) is as 

dx- f 



e{-tV ■ L°)f(X, 0) = /(X - *, 0) 



(A4) 



This simplification in Eq. (|16D yields 



/■oo 


dX~ T 


/ dr 


Jo 


dX 



((V-L^X^V-r-L^-^t-T))) 



dX 



dX- 



P . 



(A5) 



Now to express the Jacobian, X~ T and V_ r in terms of V and X we solve Eq . (|A2|) for 
short time (this is consistent with the assumption that the fluctuations are rapid [f25fl ). 
We now write the solution of Eq . (|A2j) [using Eqs.(4-6)] as follows ; 



\X N T j 



( 



X 



\ 



AT+1 



X- 



2N 



+ 



yX NJ 



( 



GtiX) 
G N (X) 



(A6) 



and 



/ y-T \ 



o1 r 



1 X ^ 



y X2N j 



( 



\ 



Gn+\{X) 



G2n{X) 



\ 



\ 



G2n{X) 



\ 



(A7) 



Here the terms of 0(t 2 ) are neglected. Since the vector X T is expressible as a function of 
X we write 



X- T = G(X) 



and the following simplification holds good; 



L\X- T , t - r).V_ T = L\G{X),t - r).V_ T 



5>l(G(X),t 



d 



T 



dXZ T 



EE L K G ( X )^-r)g 



d 



j k 



OX, 



where 



9jk 



9X, 
dXZ 



In view of Eqs.(|A6[) and ( |A7| ) we note: 



(A8) 



■j,k = l---2N (A9) 



if j — k then g^ = 1, k = 1 ■ ■ ■ N 

= e - 1T , k = N +1---2N 

if j 7^ k then oc — re _7T or 



Thus gjk is a function of r only. Let 



J^ = £4(G(X),f-T)<&* 



From Eqs. 



and (IA8I) we write 



L^A^V-r) = L}(G(X),*-r) =0 for« = l---iV 



(A10) 



(All) 



(A12) 
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So the conditions ([Ml] ), ( \Kl2}) and flpD imply that 



^•(X, i - t) = ^(Xx • - -X^, t - t) for j = 1 • • - N 
Rj(X, t-r) = ■ ■ ■ X 2Nl t-r) for j = N + 1 • • • 2N 



(A13) 



We next carry out the following simplifications of a 2 -term in Eq. 
relation (|l(J) to obtain 



We make use of the 



3 

d 2 



1,3 



3 0X i dX j 



where 







R- = E L K X ^q^ R 3 



The conditions ( |A12 ) and ( |A13| ) imply that 



(A14) 



(A15) 



R'j = for j = 1 ■ ■ ■ N 

R'j = Bl^Xi ■ • -X N ,t - t) ^ 0; for j = X + 1 ■ ■ ■ 2N 



(A16) 



By ( |A16| ) one has 



R'.VP(X,t) = V.R'P(X,t) 



(A17) 



Making use of Eqs.flTUD, (EH), fl^H ) and (TAT71) in Eq.(|A|) we obtain the Fokker-Planck 
equation fll7|). 
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FIGURES 

FIG. 1. The first four dimentionless cumulants Ai,A2,As and A4 are plotted against dimen- 
tionless time for the dynamical system described by Eq.(39). 

FIG. 2. The diffusion coefficients calculated numerically (marked as dark squares) using Eqs. 
(38) and (39) after transformation (22) are compared with theoretically obtained values (marked 
as circles) using Eq.(51) for several values of the coupling-cum-external field strength e (units are 
arbitrary) . 
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